Sub-ohmic two-level system representation of the Kondo effect 
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It has been recently shown that the particle-hole symmetric Anderson impurity model can be 
mapped onto a Z2 slave-spin theory without any need of additional constraints. Here we prove by 
means of Numerical Renormalization Group that the slave-spin behaves in this model like a two- 
level system coupled to a sub-ohmic dissipative environment. It follows that the Z2 symmetry gets 
spontaneously broken at zero temperature, which we find can be identified with the on-set of Kondo 
coherence, being the Kondo temperature proportional to the square of the order parameter. Since the 
model is numerically solvable, the results are very enlightening on the role of quantum fluctuations 
beyond mean field in the context of slave-boson approaches to correlated electron models, an issue 
that has been attracting interest since the 80's. Finally, our results suggest as a by-product that 
the paramagnetic metal phase of the Hubbard model at half-filling, in infinite coordination lattices 
and at zero temperature, as described for instance by Dynamical Mean Field Theory, corresponds 
to a slave-spin theory with a spontaneous breakdown of a local Z2 gauge symmetry. 
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Mott's localization, and all its annexes like the lo- 
cal moment formation and the Kondo effect, is a phe- 
nomenon that escapes any mean-field single-particle de- 
scription since it directly affects only part of the elec- 
trons' degrees of freedom, namely their charge. This 
is unfortunate because mean-field theory is the simplest 
and straightest way to approach interacting many-body 
systems. A trick to circumvent this difficulty, which has 
provided lots of physical insights along the years, is to 
artificially enlarge the Hilbert space adding new degrees 
of freedom that aim to describe just the electron charge 
configurations, supplemented by local constraints that 
project the enlarged Hilbert space onto the physical one. 
The final scope is to make Mott's localization accessi- 
ble already at the mean-field level. The most famous 
realization of this work-programme is the so-called slave 
boson technique, originally introduced to describe Ander- 
son and Kondo models for /-electron systems*^ 2 - A deli- 
cate issue of the slave-boson theory is that the mean-field 
treatment, although providing quite satisfactory results, 
explicitly breaks a local U(l) gauge symmetry, which 
cannot be broken hence requires going beyond mean-field 
to be restored^— More recently, novel approaches have 
been proposed in the attempt of reducing the dimen- 
sion of the enlarged Hilbert space, hence the redundancy 
of the representation, still maintaining the nice feature 
of making Mott localization accessible at the mean-field 
leveli^— For instance, in Ref. |8| it has been shown that it 
is sufficient to introduce additional slave Ising variables 
on each site, namely just two level systems rather than 
infinite level ones as in the original slave-boson technique, 
to account at the mean-field level for a Mott transition 
in the half-filled Hubbard model. In this new represen- 
tation the continuous U(l) local gauge symmetry of the 
slave-boson theory is replaced by a discrete Z 2 one. Alike 
in the slave-boson theory, the Z 2 slave-spin formulation 



must be supplemented by a constraint that projects the 
enlarged Hilbert space onto the physical one. Remark- 
ably, it has been shown in Ref. ||| that the constraint is 
unnecessary in the case of an Anderson impurity model at 
particle-hole symmetry, and, equivalently, of a Hubbard 
model at half-filling in the limit of infinite coordination 
lattices, which can be mapped^ onto an impurity model 
self-consistently coupled to a bath. Specifically, given the 
Anderson impurity Hamiltonian 
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where c] creates a conduction electron while dt an im- 
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purity one, with rid = ^2 a d\d al and the Ising plus elec- 
tron model 

n z . 2 = Hbath + o x n hyh + ^(i-o- z ), (2) 

where er a , a = x,y,z, are Pauli matrices, it follows that 
at particle-hole symmetry the following identity holds: 9 
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The Ising operator o z in |(3J) can be identified with the 
electron operator 1 — 2 (rid — 1) j which has value +1 if 
rid = 1 and -1 if rid = 0, 2, hence it describes charge 
fluctuations. Furthermore, the mixed operator o x d a ac- 
tually represents the physical impurity annihilation op- 
erator. At zero temperature, the two Hamiltonians ([1]) 
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FIG. 1: The average value (a x ) as a function of the external 
symmetry breaking field h for several values of V and U = 0.1. 

and $2$ must therefore have the same ground state en- 
ergy. It can be readily shown that a mean-field fac- 
torized wavefunction | *&) =| Ising)x | electrons) for 
the model @ allows to reproduce all mean-field results 
of the slave-boson mean-field approach to the Anderson 
impurity Hamiltonian (TTJ). In fact, if we assume that 
(Ising | a z | Ising) = cos 9 and (Ising | o x | Ising) = sin6>, 
then the electronic wavefunction must be the ground 
state of the resonant level Hamiltonian 

n* = Hbath + sin 9 n hyh + — (1 - cos 9) , (4) 

with energy E(9) and hybridization width T„ = sin 2 9T, 
lower than its bare value T. Minimization of E(9) leads 
to the same result as obtained by slave-boson mean-field 
theory. 1 ' 1 For instance, if we assume that T(e) = T for 
e 6 [-D, D] and zero otherwise, assumption that we shall 
make hereafter taking D = 1 our unit of energy, then, for 
U ^> r, we find that the known mean- field result 

sm 2 9 = ^ exp(-^ . (5) 

By analogy with slave bosons, it is tempting to interpret 
the finite value of sin 9 as manifestation of Kondo coher- 
ence, and the impurity operator d a in ([2]) as the coherent 
quasiparticle with weight sin 2 9. In fact, within mean 
field the spectral function A da (e) of the physical electron 
a x d a is simply the convolution of the spectral functions 
of the resonant level, A d (e) and of the Ising operator a x , 

A a (e) = sin 2 9 5(e) + (s(e - SI) + 5(e + fi)) , (6) 

where SI = U/2 cos 9 is the effective Zeeman splitting of 
the Ising spin. Specifically, one finds that 

A da {e) (A d * A^j (e) = sin 2 9A d {e) (7) 
+ —^- (e(e - SI) A d (e - O) + 9(-e - SI) A d (e + SI) J , 
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FIG. 2: Square of the symmetry breaking order parameter 
and of Tk/T plot in logarithmic scale versus U/T. Note the 
linear behavior at large U /Y. 

The mean-field expression of A da (e) displays a low energy 
resonance with width T* and weight sin 2 9, the rest of 
its weight being concentrated in two symmetric peaks at 
energies ±Sl. It is quite reasonable to regard these latter 
as the Hubbard side-bands and instead the central peak 
as the Abrikosov-Suhl resonance, hence T* the Kondo 
temperature Tk- However, while the mean field result 
predicts that the Hubbard bands are broadened on the 
same scale T* as the central resonance, in reality their 
width is controlled by the bare T. 

Another questionable aspect of the mean field solution 
with (a x ) ^ is that it explicitly breaks the discrete 
symmetry of ©: d a — >• —d a and a x — > —a x . Therefore, 
just like in the conventional slave-boson approach^ we 
may wonder how reliable is mean field. We are going 
to show that, unlike in slave-boson theory, the breaking 
of the discrete Z-i symmetry of © is not an artifact of 
mean-field but does spontaneously occur in the actual 
ground state. 

We start noticing that the model ^ resembles a two- 
level system coupled to a dissipative bathji 2 - where the 
two levels are the states with a x — ±1 and the tunneling 
is provided by a z . In reality, the two models are not rig- 
orously equivalent, since the hybridization operator Hhyb 
does not behave exactly like the coupling to a dissipative 
bath of bosons. However, if we judge solely from the low 
frequency behavior of the dissipative-bath spectral func- 
tion, J(uj) ~w s ^we should conclude that the model ([2]) 
behaves like a spin-boson Hamiltonian with an extremely 
sub-ohmic dissipation, s = 0. In such a case, we expect 
the tunneling to be irrelevant at low energy, hence unable 
to split the degeneracy of the two levels. This conclusion 
is remarkable because it means that the mean field re- 
sult is more correct than we would have expected; the Zi 
symmetry is spontaneously broken at zero temperature. 
There is still another aspect that we find worth mention- 
ing. In reality, the Anderson impurity model ([1} con be 
mapped into another spin-boson models with the role 
of the spin being played now by the physical spin of the 
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FIG. 3: Spectral function of the physical impurity electron 
operator, a x d a , at U = 0.1 and T = 0.01. The dots are the 
values of the impurity spectral function as calculated directly 
by the Anderson impurity model Q. We just plot few of these 
points, since the two spectral functions are just coincident. 



FIG. 4: Left panel: Spectral function of the impurity operator 
d a . Right panel: Spectral function of the Ising operator o x , on 
a large scale, top, and zoomed close to zero energy, bottom. 
The parameter that are used are U = 0.1 and V = 0.01 in 
units of the conduction electron half-bandwidth. 



impurity, rather than by the charge as in ([2]). In this 
alternative and more familiar representation, the bath is 
however ohmic and leads to incoherent derealization of 
the physical spin, once again the Kondo effect. We find 
quite amusing that the Kondo effect, known to occur in 
the Anderson impurity model ([T]) with constant T, hap- 
pens to be described by two complementary spin-boson 
models, in one of which it corresponds to derealization 
while in the other, the model Eq. ((2J, to localization. 

The above speculation can be supported by actual cal- 
culations, affordable because of the reduced dimension 
of the Hilbert space as opposed to conventional slave 
bosons. In particular, we shall consider both the An- 
derson impurity model, Eq. (JTJ, and the model and 
study them by means of the Numerical Renormalization 
Group (NRG)J^ 

The first task we undertake is to confirm that sponta- 
neous breaking of Z2 does occur in model (|5J). We hence 
add to H.Z2 a symmetry breaking perturbation —ha x , 
and calculate the average (a x ) as function of h. The 
results are shown in Fig. [TJ where it is evident the char- 
acteristic behavior of a symmetry broken phase. The 
asymptotic value of (a x ) as h — > + can be considered as 
the zero-field order parameter. In Fig. [2] we plot (a x ) 2 
as well the actual Kondo temperature in units of Y of 
the Anderson impurity model P on a logarithmic scale 
as function of U /Y. We showed that within mean field 
these two quantities coincide, but in reality they do not, 
see Fig. [2j even though the linear slope at large U/Y 
is the same, and twice as bigger as the mean field value 
— wU/16Y, Eq. (|5|). This simply means that mean field as 
usual overestimates the value of the symmetry breaking 
order parameter, hence of Tk- 

The next important step is showing that gauge invari- 
ant quantities in both models ([I]) and ((2|) are indeed the 
same. We then calculate the impurity spectral function 
of the Anderson impurity model and compare it with the 



spectral function of the physical electron a x d a in model 
@. The two quantities are found to be practically indis- 
tinguishable from each other, see Fig. [31 thus confirming 
the validity of the mapping. 

We previously showed that within mean field theory, 
the spectral function of the physical electron Ad a (e) is ap- 
proximated by the convolution between the spectral func- 
tions A a (e) of a x , whose mean field expression is given in 
Eq. ([6]), and Ad(e) of d al in mean field simply the spec- 
tral function of a resonant level model of width sin 2 9Y. 
The actual NRG Ad(e) is found not to differ much from 
the mean field result; it is still made up of a single reso- 
nance at the chemical potential, see Fig. [4[ of width the 
true Kondo temperature. Instead, A a (e) is substantially 
different from mean field, see Fig. 0] It still displays a 
(5-peak at e = with weight (cr x ) 2 , even though not ex- 
actly coincident with the value extracted in the zero-field 
limit, see the comment below. However, the finite energy 
peaks shift at much higher energy, around the edge of 
the particle-hole continuum, and get quite broadened. 
In addition, a tiny linear in e component appears at low 
energy, which resembles much the particle-hole spectrum 
of the resonant level since the linear behavior stops just 
around the Kondo temperature. This result shows that 
quantum fluctuations couple strongly the Ising spin with 
the electrons, providing a very short lifetime to the Ising 
spin-flip excitations. Although NRG, at least in the ver- 
sion we use, is not expected to be very accurate at high 
energy, especially in such a case of a sub-ohmic bath^ 
we believe that the gross features of A a (e) that we find, 
e.g. the very broad incoherent background peaked at 
high energy, are true. 

We observe that the spectral function A<fo-(e) of the 
physical electron o x d G1 shown in Fig. [31 can be generally 
written as 

A da (e) = fdu} A d (e - u) A a (cj) K{u, e), (8) 
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FIG. 5: Convolution of the spectral functions A a {e) and Ad(e) 
for U = 0.1 and T = 0.01. Note the absence of Hubbard 
bands, which should be appear at ~ ±(7 /2 = 0.05 but are hid- 
den below the broad background peaked around the particle- 
hole band edge at energy ~ ±2D = ±2. 

where the kernel K(u), e) amounts for all vertex correc- 
tions. In Fig. [S]we show the simple convolution of Ad(e) 
and A a {ui), which, compared with the correct result in 
Fig. [3l could provide a rough estimate of vertex correc- 
tions. We notice that, while the Abrikosov-Suhl reso- 
nance is reproduced quite well also by the convolution, 
the Hubbard side-bands are completely masked by the 
broad background of A a (e), see Fig. [4] This implies that 



vertex corrections play a major role at high energy, fil- 
tering out that high-energy background and letting the 
Hubbard bands emerge. 

In conclusion, we have shown several amusing fea- 
tures of the i?2 slave-spin representation of the particle- 
hole symmetric Anderson impurity model that, unlike its 
slave-boson analogous, can be exactly solved by NRG. 
Specifically, we have shown that in this language the 
Kondo effect corresponds to spontaneous breaking of a 
local discrete Z2 symmetry, that does survive quantum 
fluctuations. 

We end by mentioning that, as a byproduct, our result 
suggests that the zero-temperature paramagnetic metal 
phase of the half-filled Hubbard model, in the limit of in- 
finite lattice-coordination^ can be regarded within the 
Z2 slave-spin representation 8 '- as a phase where a local 
Z2 gauge symmetry is spontaneously broken. Conversely, 
the zero-temperature Mott metal-to-insulator transition 
would correspond to the restoration of the Z2 symme- 
try. We note that this result does not violate Elitzur's 
theorem,— because of the infinite lattice-coordination 
limit »i£ However, we expect that the local Z2 gauge sym- 
metry should be recovered at any finite temperature, 
since a sub-ohmic dissipative two-level system is known 
to delocalize at any finite temperature^ 
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